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Objective Metabolic abnormalities, e.g., diabetes, are common among schizophrenia patients. Peroxisome proliferator activated 
receptor-Y (PPAR-y) regulates glucose/lipid metabolisms, and schizophrenia like syndrome may be induced by actions involving reti- 
noid X receptor-a/PPAR-Y heterodimers. We examined a possible role of the PPAR-y gene in metabolic traits and psychosis profile in 
schizophrenia patients exposed to antipsychotics. 

Methods Single nucleotide polymorphisms (SNPs) of the PPAR-y gene and a serial of metabolic traits were determined in 394 schizo- 
phrenia patients, among which 372 were rated with Positive and Negative Syndrome Scale (PANSS). 

Results SNP-10, -12, -18, -19, -20 and -26 were associated with glycated hemoglobin (HbAlc) whereas SNP-18, -19, -20 and -26 were 
associated with fasting plasma glucose (PPG). While SNP-23 was associated with triglycerides, no associations were identified between 
the other SNPs and lipids. Further haplotype analysis demonstrated an association between the PPAR-y gene and psychosis profile. 

Conclusion Our study suggests a role of the PPAR-y gene in altered glucose levels and psychosis profile in schizophrenia patients ex- 
posed to antipsychotics. Although the Prol2Ala at exon B has been concerned an essential variant in the development of obesity, the 
lack of association of the variant with metabolic traits in this study should not be treated as impossibility or a proof of error because oth- 
er factors, e.g., genes regulated by PPAR-y, may have complicated the development of metabolic abnormalities. Whether the PPAR-y 
gene modifies the risk of metabolic abnormalities or psychosis, or causes metabolic abnormalities that lead to psychosis, remains to be 
examined. Psychiatry Investig 2014;11(2):179-185 
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INTRODUCTION schizophrenia patients. DM may be induced by antipsychotics 

medications,' or may occur before the antipsychotics exposure 
Metabolic abnormalities, e.g., DM, are common among in a considerable number of schizophrenia patients.^ Positive 
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syndrome of most schizophrenia patients can be relieved to 
some extents by antipsychotic treatments, but this is not the 
case for negative syndrome,^''' suggesting that schizophrenia 
patients with severe negative syndrome may represent a sub- 
group with distinct biological features. 

Peroxisome proUferator- activated receptors (PPARs) belong 
to the superfamily of nuclear receptors (NRs), and play a cen- 
tral role in the transcriptional control of glucose and lipid reg- 
ulation. The agonism of NRs has been used as a molecular 
strategy for treating diabetes and hyperglycemia."^ Three dif- 
ferent human PPAR subtypes, namely PPAR-a, PPAR-(3 (also 
known as PPAR-6) and PPAR-Y,*"'' have been identified. 
PPAR-y is the most frequently studied NR, and is involved in 
the control of energy balance and Upid and glucose homeo- 
stasis.* Several PPAR-y gene knockout studies in mice have 
demonstrated glucose and lipid metabolism cross-talks be- 
tween adipose tissue, muscle and Uver, suggesting that PPAR-y 
is important for retaining normal insulin sensitivity and 
whole-body glucose and lipid homeostasis.'''" Also, genome- 
wide association and association studies have indicated a ge- 
netic contribution of the PPAR-y gene in the development of 
type 2 DM and metabolic traits." '^ 

The activation of PPAR-y mediates anti-inflammation in 
the cardiovascular system and nervous system." PPAR-y 
forms a heterodimer complex with retinoid X receptor-a 
(RXRa), which then binds to peroxisome proliferator re- 
sponse elements (PPREs) within the promoters of PPAR-y- 
targeted genes. " ''^ Retinoic acid receptors (RARs) are present 
in all parts of the cranial region, and delivery of retinoids is 
exquisitely controlled throughout the embryonic and fetal de- 
velopment."" RARs act either alone or together with the het- 
erodimeric partners. An in vitro experiment in adult human 
keratinocytes showed that endogenous retinoic acid RAR- 
RXR heterodimers are the major functional forms that regu- 
late the retinoid-responsive elements,'' suggesting that an ac- 
curate regulation between PPAR-y, RAR and RXR is needed 
to ensure a normal brain development and function. Because 
retinoid toxicity or deficit produces symptoms resembling 
schizophrenia,"' it is reasonable to deduce that genetic variants 
in the PPAR genes may affect the functioning of retinoid cas- 
cade and, thus, have profound implications in schizophrenia. 

The objective of the study was to investigate a possible role 
of the PPAR-y gene in metabolic traits and psychosis profile 
in schizophrenia patients in Taiwan. 

METHODS 
Participants 

This study was approved by the Institutional Review Board 
of Yu-Li Hospital, Department of Health, Taiwan, file number 



YLH-93001. A total of 401 subjects (300 males and 101 fe- 
males) meeting DSMIV criteria of schizophrenia and who 
had been exposed to antipsychotics treatment for more than 
2 years, were recruited. None of the patients were exposed to 
psychiatric therapeutics other than antipsychotics except hyp- 
notics when medically required. Medication for physical 
health problems was intervened in any form during the study 
period. The mean age of the study population was 49.3+9.9, 
ranging from 21 to 84 years. AU eligible patients signed the 
informed consent form for blood donation and genetic analy- 
sis. At last, 394 good satisfactory DNA samples were obtained 
and used for genotyping. Psychosis profiles were evaluated 
using the Positive and Negative Syndrome Scale (PANSS) Tai- 
wan Chinese version, which includes subscales of positive 
syndrome, negative syndrome and general psychopathology.'* 
Three hundred and seventy two out of the 401 schizophrenia 
patients were successfully assessed using PANSS by board- 
certified psychiatrists. None of the patients were rated during 
acute psychosis. All patients examined in this study resided in 
a commodious environment. 

Metabolic traits 

BMI, waistline, fasting plasma glucose (FPG), glycated he- 
moglobin (HbAic), total cholesterol (TC), triglycerides (TGA), 
low density lipoprotein cholesterol (LDL), high density lipo- 
protein cholesterol (HDL), were measured on October 2004. 
Details of measurement of the indexes are summarized in 
supplement file. 

Genotyping 

The human PPAR-y gene is located on chromosomal re- 
gion 3p25." At least 7 human PPAR-y mRNA splice variants 
have been found, each consists of exons 1 through 6 consecu- 
tively.^" The PPAR-y protein exists in 2 isoforms-PPAR-yl 
with 477 amino acids and PPAR-Y2 with 505 amino acids. 
The PPAR-yl mRNAs are coded by 8 exons (Al, A2, 1, 2, 3, 4, 
5, and 6), among which exon Al and A2 are untranslated. The 
PPAR-Y2 mRNAs are encoded by 7 exons (exon B, 1, 2, 3, 4, 5, 
and 6), among which the additional 28 amino acids are en- 
coded by exon B.^' 

Information of single nucleotide polymorphisms (SNPs) of 
the PPAR-y gene in Han Chinese population (Beijing, China) 
and Japanese population (Tokyo, Japan) were retrieved from 
the HapMap database (phase II, 7 Jan 2007, http://hapmap. 
ncbi.nlm.nih.gov). Twenty TaqSNPs with minor allele fre- 
quency (MAP) >0.03 were selected using the Hapmap Tagger 
Program; 9 coding SNPs available during the research period 
were selected from dbSNP (Build 125) for genotyping. De- 
tailed information of the SNPs, including chromosomal posi- 
tion, nucleotide location, exonic or intronic (based on PPAR- 
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yl mRNA and protein encoded by exon B, 1, 2, 3, 4, 5, and 6), 
allelic variants and amino acid changes, are summarized in 
supplement file. The Sequenom MassARRAY* iPLEX Gold 
platform (Sequenom Inc., San Diego, CA, USA) and caUing 
software Typer were used for genotyping and signal analysis. 

Statistical analysis 

SNPs not deviated from the Hardy- Weinberg equilibrium 
and >0.03 were analyzed using the Haploview™ program 4.1 
(http://www.broadinstitute.org/ scientific-community/ sci- 
ence/programs/medical-and-population-genetics/haploview/ 
haploview). General linear model (GLM) was used to esti- 
mate the contribution of each SNP, age, gender, the use of 
typical antipsychotics, and the use of atypical antipsychotics 
on metabolic indexes by applying the enter method. SNPs 
that possibly confer fiinctionality or associated with genetic 
variants with fiinctionality (i.e., SNPs significantly associated 
with metaboUc indexes) were used to evaluate a role of the 
PPAR-y gene in psychosis profile. The additive and dominant 
genetic efiects of each haplotype on PANSS were tested against 
the other haplotypes (i.e., other haplotypes were merged and 
considered as one haplotype) in the same GLM with age, gen- 
der, the use of typical antipsychotics, and the use of atypical 
antipsychotics included. The code for additive effect (0, 1 or 2) 
in the data sheet was the frequency of the haplotype of inter- 
ests. For dominant effects, the code was 1 when one haplotype 
of interests existed, and 0 when no haplotype of interest exist- 
ed. After each haplotype was examined, those significantly 
contributed to PANSS either through additive or dominant ef- 
fects, were further permutated with 10,000 iterations. The per- 



mutation tests randomly assigned the identification labels to 
the patients (each possessed a specific haplotype), creating a 
scenario where the dependent variable (i.e., PANSS) was not 
related to the PPAR-y gene. The SAS program was used for the 
GLM statistical analysis and permutations. 

RESULTS 

Among the 29 SNPs we examined, genotype distribution of 
2 SNPs significantly deviated from the Hardy-Weinberg equi- 
librium and 9 SNPs had MAP <0.03 (Supplementary Table 1, 
online). These 1 1 SNPs were excluded from statistical analysis 
to avoid a bias toward significance. Figure 1 shows the linkage 
disequilibrium (LD) strengths in D' between the 18 remained 
SNPs and haploblocks revealed by the Haploview program. 
The LD plot consisted 3 haploblocks: block- 1 included SNP-1 
(rs6785890, 5'UTR) and -2 (rsl2631819, 5'UTR), spanning a 
7 kb region; block 2 included SNP-4 (rsl0510410, 5'UTR), -5 
(rsl3061415, 5'UTR), -6 (rsl2485478, 5'UTR), -9 
(rsl2496005, 5'UTR), -10 (rsl2490265, 5'UTR) and -11 
(rsl801282, Prol2Ala, exon B), spanning a 46-kb region, fol- 
lowed by the unblocked SNP- 12 (rsl7817276, intron B), -14 
(rsl3306745, intron B) and -17 (rs4135268, intron B); block-3 
included SNP-18 (rs2959272, intron 3), -19 (rsl875796, in- 
tron 3), -20 (rs4135275, intron 3), -21 (rsl2489347, intron 3), 
-23 (rsl3306747, Pro269Pro, exon 5) and -26 (rsl 175544, in- 
tron 5), spanning a 24 kb region, followed by the independent 
SNP-29 (rs3856806, His449His, exon6). None of the SNPs 
were located at exon 1, intron 1, exon 2, intron 2 or exon 3. 
Further information on the haplotype frequencies and 
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Figure 1. LD between valid SNPs. LD: 
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strengths between blocks are shown in Figure 2. The LD be- 
tween block- 1 and block-2 was 0.89, whereas that between 
block-2 and block-3 was 0.82. The unblocked SNP-12, -14 
and -17 were located in intron B, suggesting that numerous 
historical recombinations occurred between intron B and 
block-2 and block-3. In terms of molecular evolution, block- 1, 
block-2 and block-3 were possibly derived from the same 
DNA segment, while intron B was derived from DNA else- 
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Figure 2. LD between haploblocks. LD: Linkage disequilibrium. 



where. 

Significance of SNPs by allele distribution without consid- 
ering gender interaction is summarized in Figure 3. SNP-10, 
-12, -18, -19, -20, and -26 were significantly associated with 
HbAic, whereas SNP-18, -19, -20, and -26 were associated 
with FPG. Only SNP-23 was significantly associated with 
TGA. No associations were found between the SNPs and oth- 
er metabolic indexes under the same parameters. Significance 
of SNPs by allele distribution when gender interaction was 
included is summarized in Figure 4. Similar associations were 
observed for HbAic, FPG and TGA. 

Six SNPs significantly associated with HbAic including 
SNP-10, -12, -18, -19, -20, and -26, were selected to form hap- 
lotypes and used in the GLM and permutation tests. SNP-23 
was tested separately because it was urrrelated to glucose regu- 
lation in our statistical model. Eleven haplotypes were identi- 
fied using the SAS program (Supplementary Table 2, online). 
Each haplotype was regressed with PANSS adjusted for age, 
gender, the use of typical antipsychotics and the use of atypi- 




Figure 3. Significance of the SNPs at the 
PPAR-y gene in metabolic indexes evalu- 
ated by allele distribution, without gender 
interaction. P-values are shown in the 
form of minus natural logarithm (-In P). 
The green line is the significant level of 
-In 0.05 (x^=2.99). Rs serial numbers cor- 
responding to the SNPs are summarized 
in Supplementary Table 1. SNPs: single 
nucleotide polymorphisms, PPAR-y: per- 
oxisome proliferator activated receptor-y. 
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Figure 4. Significance of the SNPs at the 
PPAR-y gene in metabolic indexes eval- 
uated by allele distribution, with gender 
interaction. SNPs: single nucleotide poly- 
morphisms, PPAR-y: peroxisome prolif- 
erator activated receptor-y, BMI: Body 
Mass Index, FPG: fasting plasma glu- 
cose, HbAic: glycated hemoglobin, TC 
total cholesterol, TGA: triglycerides, LDL 
low density lipoprotein cholesterol, HDL 
I high density lipoprotein cholesterol. 
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Table 1. Statistical results for GLMs and permutation tests (10,000 iterations) 



Haplotype 


Alleles 


Frequency 




GLM 






Permutation test 




T 


PS 


NS 


GP 


T 


PS NS 


GP 


6 


G-A-G-T-G-T 


0.0095 


















Additive 




0.3821 


0.0064* 


0.2615 


0.1790 




0.0284* 






dominant 


















7 


G-A-T-C-A-C 


0.0157 


















Additive 




0.0015* 


0.0058* 


0.1259 


0.0007* 


0.0091* 


0.0192* 


0.0068* 




dominant 


















8 


G-A-G-T-A-C 


0.1005 


















Additive 




0.0134* 


0.0012* 


0.4219 


0.0156* 


0.0446* 


0.0138* 


0.0505 




dominant 




0.1004 


0.0122* 


0.4952 


0.2048 




0.0409* 





*significant p- values. T: total scale, PS: positive syndrome, NS: negative syndrome, GP: general psychopathology, GLM: General linear model 



cal antipsychotics in the GLM. The results are summarized in 
Table 1. Haplotype 6 and 7 were significantly associated with 
one or more scales in additive mode only. Haplotype 6 was 
associated with positive syndrome (p=0.0064) but not with 
total scale (p=0.3821), negative syndrome (p=0.2615) or gen- 
eral psychopathology (p=0.1790). Haplotype 7 was associat- 
ed with total scale (p=0.0015), positive syndrome (p=0.0058), 
and general psychopathology (p=0.0007), but not with nega- 
tive syndrome (p=0.1259). Haplotype 8 was associated with 
total scale (p=0.0134), positive (p=0.0012) and general psy- 
chopathology (p=0.0156) in additive mode, but no with nega- 
tive syndrome (0.4219); the haplotype also demonstrated an 
association with positive syndrome (0.0122) in dominant 
mode. These 3 haplotypes, with a summed up frequency of 
12.6%, were then examined by permutation tests. All associa- 
tions remained significant after 10,000 iterations, except for 
the association of haplotype 8 with general psychopathology 
(p=0.0505). The same alleles, namely allele G and A, existed 
in SNP-10 and 12 in the 3 haplotypes. 

The regression of SNP-23 in PANSS demonstrated no evi- 
dence of association (data not shown). 

DISCUSSION 

SNP-ll, the Prol2Ala nonsynonymous variant at exon B, is 
one of the most widely examined SNP at the PPAR-y gene in 
studies investigating obesity related phenotypes, including in- 
sulin resistance and type 2 diabetes."''^^ When this SNP was 
used to interpret the results, PPAR-Y2 mRNA or protein were 
presumed to play a major role over PPAR-yl mRNA or pro- 
tein in the phenotype of interests, because exon B exists only 
in the PPAR-Y2 isoform. PPAR-yl mRNA is predominant in 
human tissues, but PPAR-Y2 mRNA only appears in human 
Uver and adipose tissues. We found no evidence supporting a 



major role of the Prol2Ala at exon B (and thus PPAR-Y2) in 
maintaining a normal glucose level. Although the Prol2Ala 
at exon B has been concerned an essential variant in the de- 
velopment of obesity, that lack of association of the variant 
with metabolic traits in this study should not be treated as 
impossibility or a proof of error because other factors, e.g., 
genes regulated by PPAR-y, may have complicated the devel- 
opment of metabolic abnormalities. SNP-23 (Pro269Pro), a 
synonymous variant at exon 5, was associated with TGA by 
allele distribution. Although a role of the PPAR-y gene in lipid 
profile in schizophrenia patients remains inclusive, fiirther in- 
vestigation is worth performing, especially the lipid pathway 
concerning the Apolipoprotein E (APOE). The APOE gene is 
a schizophrenia gene.^^'^*" This gene is also a diabetic gene, 
which confers its risk at least partially through altering serum 
lipid levels.^' Importantly, APOE is transcriptionally regulated 
by PPAR-y and liver X receptors (LXRs),^* which form obligate 
heterodimers with RXRs. Variations in the PPAR-y gene may 
change the risk schizophrenia, psychosis syndrome or diabe- 
tes, through the subtle actions of PPAR-y:RXR and LXR:RXR. 
For the haplotype effects on the PANSS scores, haplotype 6, 7 
and 8 demonstrated a positive effect on positive syndrome al- 
though their frequencies were low. The common GA alleles of 
the two SNPs embracing exon B may, or are in linkage with 
genetic variants that affect the severity and long-term out- 
come of schizophrenia. 

There are a few limitations in this study. First, our study 
subjects were schizophrenia patients undergoing long term 
antipsychotic treatment. The proclivity of antipsychotic drugs 
for metabolic disturbances, especially the second-generation 
(atypical) antipsychotics, have been extensively documented. 
For this reason, the associations identified in this study might 
be caused by the adverse effects of antipsychotics rather than 
the disease per se. Nevertheless, contributions of the antipsy- 
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chotics were accounted for in the regression analysis. Third, 
haplotype 6 and 7 were of very low frequencies (<1% and 
1.6%, respectively) and thus associations of these two genetic 
structures with PANSS relied on very small number of pa- 
tients. Nevertheless, haplotype 8 was of relative high frequen- 
cy (10.1%), suggesting the association was likely a true associ- 
ation. 

Schizophrenia is a heterogeneous disorder in etiology, symp- 
tomology and treatment responses.^' Albeit metabolic abnor- 
malities are associated with cognitive dysfiinctions and brain 
abnormalities, so far, there have been limited evidences show- 
ing that metabolic abnormalities lead to specific mental disor- 
ders.'" The present study demonstrates a significant associa- 
tion of the PPAR-y gene in altered glucose levels and psychosis 
profile in schizophrenia patients. Whether the PPAR-y gene 
modifies the risk of metabolic abnormalities and psychosis, or 
causes metabolic abnormalities that lead to psychosis, remind 
to be determined. 
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